A phason disordered two dimensional quantum antiferromagnet 
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We examine a novel type of disorder in quantum antiferromagnets. Our model consists of localized 
spins with antiferromagnetic exchanges on a bipartite quasiperiodic structure, which is geometrically 
disordered in such a way that no frustration is introduced. In the limit of zero disorder, the structure 
is the perfect Penrose rhombus tiling. This tiling is progressively disordered by augmenting the 
number of random "phason flips" or local tile-reshuffling operations. The ground state remains Neel 
ordered, and we have studied its properties as a function of increasing disorder using linear spin 
wave theory and quantum Monte Carlo. We find that the ground state energy decreases, indicating 
enhanced quantum fluctuations with increasing disorder. The magnon spectrum is progressively 
smoothed, and the effective spin wave velocity of low energy magnons increases with disorder. For 
large disorder, the ground state energy as well as the average staggered magnetization tend towards 
limiting values characteristic of this type of randomized tilings. 

PACS numbers: 71.23.Ft, 75.10.Jm, 75.10.-b 



This paper discusses the effects of a novel type of ge- 
ometric disorder on the ground state of Heisenberg anti- 
ferromagnets. The clean system is a perfectly determin- 
istic quasiperiodic two-dimensional tiling, and the type 
of disorder we consider is geometric, involving a discrete 
shift of a randomly selected subset of sites. These pha- 
son flips are operations that correspond to reorganizing 
the structure locally, in the vicinity of the flipped site. 
This type of disorder is strongly constrained, as one does 
not modify the basic building blocks of the structure, 
but only the way they are connected. For many qua- 
sicrystalline alloys, samples of good quality are believed 
to possess long range order, with some thermally induced 
phonon and phason disorder present at finite tempera- 
tures However, one can also consider essentially ran- 
dom quasiperiodic structures, which could be preferred 
for entropic reasons Q. 

Quenched phason disorder has been considered pre- 
viously in other contexts, in particular concerning elec- 
tronic properties, and the possibility of Anderson local- 
ization in such systems was discussed. Benza et al. Q 
considered quantum diffusion in a two dimensional ran- 
domized tiling, and found that the typical value of the 
diffusion exponent was larger for the disordered system 
compared to the pure case. This means that quasicrys- 
tals show delocalization-like effects from disorder rather 
than the opposite, weak localization due to disorder, ob- 
served in periodic systems. Piechon and Jagannathan 
came to the same conclusion by analyzing the statistics 
of the energy levels in phason disordered tilings (see e.g. 
Refs. for reviews). Schwabeetal. 0] analyzed the 

effect of phason flips on the electronic levels and wave- 
functions, showing that they lead to a smoothing of the 
density of states, as well as of the fluctuations of the 
conductance. Finally, the effects of randomness on the 




FIG. 1: Illustration of a perfect (upper panel) and a phason 
disordered Penrose tiling (lower panel). 



phonon spectrum has been considered in two dimensional 
quasiperiodic tilings (see e.g. Ref. |8j). 

The study of the Heisenberg model on a quasiperiodic 
tiling is motivated by experimental findings of antifer- 
romagnetic correlations in quasiperiodic ZnMgR alloys 
(R: rare earth) [9fl . fn theoretical models considered thus 



far [lCj, |ll|, |rj, |13, bipartite tilings are considered, 



with all sites occupied by spins, and their exchanges re- 
stricted to adjacent lattice sites, such that no frustration 
arises. Furthermore, these studies did not take into ac- 
count the effects of disorder, which is almost certainly 
also present in the experimentally studied alloys. Here, 
we focus on the effects that arise in quasiperiodic mag- 
nets in the presence of phason disorder, which is found 
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FIG. 2: (Color online) A single phason flip, showing the origi- 
nal (ro) and final (r' ) position of the shifted lattice site. New 
edges are shown by dashed lines. 



to lead to enhanced quantum fluctuations, as indicated 
by the lowering of the ground state energy along with a 
reduction of the staggered magnetization with increasing 
disorder. 

Fig. Q] shows a portion of a perfect (deterministic) 
Penrose tiling and a typical example of a phason dis- 
ordered tiling. We consider in the following the antifer- 
romagnetic spin-1/2 Hcisenberg model with Hamiltonian 
H = Jj2(ij) S« • Sj for spins Si located on all vertices 
of such tilings. Nearest-neighbor exchanges are antifer- 
romagnetic J > 0, and act between pairs of sites that are 
linked by an edge. Despite the fact that all couplings are 
equal, the ground state is spatially inhomogeneous con- 
trarily to e.g. the square lattice antiferromagnet, where 
the staggered magnetization is uniform. The local stag- 
gered magnetizations vary as a function of the local en- 
vironment, and the ground state takes a complex, hier- 
archically organized structure 0, EH, H Ei El- The 

disorder considered here is purely geometric, i.e. the cou- 
pling J along bonds remains fixed at a constant value. 

Using linear spin wave theory (LSWT) and quantum 
Monte Carlo (QMC) calculations, we consider periodic 
approximants of the Penrose tiling Ell- These are 
finite samples of N spins, satisfying periodic boundary 
conditions. Perfect tilin gs a re obtained by the cut-and- 



project method (see e.g. [14|) after which they are disor- 



dered, by the following method: A phason flip is a process 
by which a 3-fold site hops to a new allowed position (in 
terms of the tile configurations). The old site disappears, 
as do the three bonds linking it to its neighbors, while a 
new 3-fold site appears on the other sublattice (Fig. [2]). 
In our phason generating procedure, we randomly select 
a three-fold site. If ro and fl (j = 1, 2, 3) denote the posi- 
tion vectors of the central site and its three neighbors, the 
new position of the site is given by r' — Tq = X^=i ( r j — r o) 
(Fig. [5J). Three new bonds appear linking the new site 
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FIG. 3: (Color online) Dependence of the ground state en- 
ergy on the disorder strength A calculated within LSWT. 
The smooth curve is a fit to an exponential decay towards 
a higher disorder value. The system size A=4414. The in- 
set shows the normalized ground state energy as a function 
of disorder A. Points represent values obtained after scaling 
LSWT and QMC data to infinite system sizes. 



to the sites at the positions r'j (j — 1,2,3). The coor- 
dination numbers of all of the seven sites involved are 
updated, and the whole procedure is repeated, with the 
constraint that there be an equal number of flipped sites 
on the A and the B sublattice, so as to preserve the con- 
dition Na = Nb, where Na (Nb) denotes the number of 
sites of the A (B) sublattice. 

For a given total number of phason flips N p h, the 
degree of disorder A is defined as the average overlap 
distance between the perfect sample and the disordered 
samples. The overlap for a given sample is defined by 
Nph £i=i Vii where 77, is 1 or 0, depending on whether 
site i is shifted or not with respect to the perfect tiling. 
The whole procedure is then repeated for approximants 
of system sizes N = 246, 644, 1686 and 4414 in LSWT 
and 246, 644 and 1686 in QMC. Averaging is carried out 
over a large number of samples: 100 samples for N = 246, 
644, 1686 and 10 samples for N = 4414 in LSWT and 
40 samples for N = 246 and 20 samples for N = 644, 
1686 in the QMC. This method of disordering does not 
change the overall number of rhombi of each kind, and 
generates samples of fixed phason strain [l5| . Due to the 
fact that we obtain a defected structure (a periodic ap- 
proximant, which is only locally equivalent to the infinite 
quasiperiodic tiling) there is an upper limit to the num- 
ber of phason flips we are able to introduce in the tilings. 
Unacceptable configurations are observed to appear after 
a number of flips larger than about 1.3A^, corresponding 
to the number of steps after which the periodic boundary 
conditions are felt by the system. 

After the randomized samples are obtained. LSWT cal- 
culations are performed as described in [12J, [l4|. Spin 
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FIG. 4: Evolution of the density of states p(E) for various 
values of the disorder strength A from LSWT on the N — 
4414 sites system. 
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FIG. 5: (Color online) Normalized averaged staggered mag- 
netization as a function of phason disorder strength A. Points 
represent values obtained from LSWT and QMC data after 
scaling to the thermodynamic limit. The inset shows the aver- 
aged spin wave velocity c as a function of the disorder strength 
A as obtained within LSWT. 



operators are transformed using the Holstein-Primakoff 
transformation to bosonic operators <ij, bj (i,j = 
l,...,N/2), corresponding the A and B sublattices re- 
spectively. The linearized Hamiltonian in the boson op- 
erators is then diagonalized numerically. Once the eigen- 
modes have been determined, one obtains the ground 
state energy, and local staggered magnetizations for each 
realization of disorder. Finally, we carry out the statis- 
tical analysis of the results, by performing disorder aver- 
aging. Within the QMC simulations, we obtain the lo- 
cal staggered magnetization m 2 s {i) — y.f— i e i e i (SfSp 
from the spin-spin correlation function [lCt Il2j . where 
ej ± 1, depending on whether lattice site i belongs to 
sublattice A or B. The QMC simulations were performed 
using the stochastic series expansion method [16| at tem- 
peratures taken low enough to obtain ground state prop- 
erties of these finite systems fl2| . 

We now describe our results on the effect of phason 
disorder. First, we consider the average ground state en- 
ergy per site Eq/N, shown as a function of disorder in 
Fig. [3] for the N = 4414 sites system. The values have 
been normalized with respect to the value obtained for 
the clean system. Also shown is a fit to an exponen- 
tial decay towards a limiting value in the strong disorder 
case, of the form E /N — + (e per f — edis)e~ aA , with 
a = 11.16, e dls = -0.6500 and e perf = -0.6429. The 
asymptotic value of the ground state energy edis rep- 
resents the average value of the ground state energy of 
maximally randomized Penrose approximants, which lies 
below the ground state energy of the perfect system. This 
indicates that the introduction of phasons tends to en- 
hance quantum fluctuations in the tilings, as compared 
to the clean case. The inset of Fig. [3] shows the ground 
state energy per site as a function of disorder, obtained 
from LSWT and QMC simulations, normalized with re- 
spect to the value obtained for the clean system. The 



results obtained from LSWT and QMC are found to be 
in good agreement. This also shows the applicability of 
the linear spin wave approximations to this randomized 
system. 

Next, we consider the density of states, defined via 
P(E) — S M ^(-^ — k-V)' an d shown in Fig. [4] We find that 
the dominant effect of phason disorder on the density of 
states is to progressively smoothen fluctuations and fill 
in gaps, as in the electronic case The low energy 
tail is quadratic, and can be fitted to obtain the aver- 
aged low-energy spin wave velocity c, that increases with 
A as shown in the inset of Fig. [5l This indicates that 
spin wave propagation is facilitated by the phason dis- 
order, in analogy with the problem of quantum diffusion 
of electrons in the tight binding model in quasiperiodic 
tilings [4] . In addition, the localized states at E = 3 dis- 
appear progressively. These states arise on closed loops 
of 3- fold sites They are hence destroyed when a 

phason flip occurs on one of the participating sites. 

Finally, we turn to discuss the evolution of the stag- 
gered magnetization upon introducing phason disorder. 
Fig. [5] shows the spatial average of the staggered mag- 
netization (i.e. after averaging over all of the sites) as a 
function of disorder, normalized with respect to the value 
obtained for the clean system. The curve shows a clear 
decrease in the global staggered magnetization with in- 
creasing disorder. As for the ground state energy curve, 
the decrease eventually levels off. It is also interesting to 
analyze the evolution of the full distribution of the stag- 
gered magnetizations with increasing phason disorder. In 
the perfect tiling, this distribution has several peaks, each 
of which corresponds to a distinct coordination number. 
In the disordered tilings, the differences between the co- 
ordination number is smoothed out. As Fig. [5] shows, 
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FIG. 6: Evolution of the distribution in the local staggered 
magnetization with the disorder strength, as obtained within 
LSWT on the 4414 sites system. 
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FIG. 7: (Color online) Color representation of the local stag- 
gered magnetizations in a perfect tiling (left) and a disor- 
dered tiling (right) showing sites of magnetizations ranging 
from high (blue) to low (red) values. 



locity in the disordered tiling as compared to the perfect 
quasiperiodic tiling. The density of states is smoothed. 
Degenerate states localized on closed loops disappear 
with increasing disorder. Eigenmodes tend to become 
more delocalized as compared to the perfect tiling. These 
effects have their analogy in phonon models, as well as 
in the tight binding model for electrons in quasiperiodic 
tilings. In the magnetic problem our results also indi- 
cate that the effect of disorder is to reduce the strong 
coherent backscattering of the magnon wavefunctions 
due to a perfect quasiperiodic potential and to favor a 
more diffusive dynamics. Upon increasing phason dis- 
order in the antiferromagnetic model, the ground state 
progressively loses its self-similar features, and tends to- 
wards a more homogeneous distribution of staggered mo- 
ments. The global average of the staggered magnetiza- 
tion decreases, as does the ground state energy, signal- 
ing increased quantum spin fluctuations. Both quantities 
tend towards a limiting value for this class of disordered 
tilings. Such two dimensional quantum antiferromagnets 
may be experimentally realizable in the near future, as 
there has been considerable progress recently with de- 
positing atoms on quasiperiodic surfaces, which serve as 
templates. A monolayer of Pb atoms having quasiperi- 
odic symmetry has been obtained [13] • To study quan- 
tum magnetism one needs to obtain such a surface layer 
using low spin atoms. In addition, one needs the inter- 
actions to be predominantly short range and antiferro- 
magnetic, and this could be realized via a superexchangc 
mechanism, as in the cuprate layers of high T c com- 
pounds, in which oxygen atoms mediate the exchange 
between the 5=1/2 copper spins. 

We thank HLRS Stuttgart and NIC Jiilich for alloca- 
tion of computing time. A.S. was supported by the Eu- 
ropean Commission, through a Marie Curie Foundation 
contract, MEST CT 2004-51-4307, during the course of 
this work. 



the distribution becomes smoother, as the disorder is in- 
creased. In addition, the average value shifts to lower 
values. The smoothing occurs due to a larger number of 
local environments created by the phason flips, and due 
to the loss of self similarity on larger length scales. 

Fig. [7] shows, that the distribution is more homoge- 
neous in the disordered tiling, which furthermore lacks 
the hierarchical features of the perfect tiling. The redder 
appearance of the right figure also illustrates the already 
noted fact, that globally the staggered magnetization is 
lower in the disordered system. 

In conclusion, we analyzed the effects of a new type 
of geometrical disorder on the magnetic properties of 
quasiperiodic antiferromagnets. Concerning the low en- 
ergy modes, magnons propagate with a slightly higher ve- 
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